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Using the spin-spiral formulation of the tight-binding linear muflin-tin orbital method, the prin- 
cipal components of the exchange stiffness tensor are calculated for typical hard magnets including 
tetragonal CoPt-type and hexagonal YC05 alloys. The exchange stiffness is strongly anisotropic 
in all studied alloys. This anisotropy makes the domain wall surface tension anisotropic. Com- 
petition between this anisotropic surface tension and magnetostatic energy controls the formation 
and dynamics of nanoscale domain structures in hard magnets. Anisotropic domain wall bending is 
described in detail from the general point of view and with application to cellular Sm-Co magnets. 
It is shown that the repulsive cell-boundary pinning mechanism in these magnets is feasible only 
due to the anisotropic exchange stiffness if suitably oriented initial pinning centers are available. 
In polytwinned CoPt-type magnets the exchange stiffness anisotropy controls the orientation of 
macrodomain wall segments. These segments may reorient both statically during microstructural 
coarsening and dynamically during the macrodomain wall splitting in external field. Reorientation 
of segments may facilitate their pinning at antiphase boundaries. 
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I. INTRODUCTION 

The formation and dynamics of magnetic domain 
structures are commonly studied by micromagnetic 
methods employing the phenomenological gradient ex- 
pansion for the 'exchange term' in the free energyi For a 
crystal of arbitrary symmetry this term may be generally 
written as 



A. 



Q/3 



de^ de.y 
dra drp 



(1) 



where e(r) = M(r)/M is the unit vector parallel to mag- 
netization, summation is assumed over repeated Carte- 
sian indices, Aap = MDap/^, and Dq/? is the spin- wave 
stiffness tensor which determines the long-wavelength 
part of the magnon spectrum as Tujj{c\) = Dapqaqp- In 
a cubic crystal Aap = A Sap where Sap is the Kronecker 
symbol, and A is commonly referred to as the exchange 
constant. Below Aap is referred to as the exchange stiff- 
ness tensor. 

A cubic crystal may only have a fourth-order magne- 
tocrystalline anisotropy (MCA) in the spin-orbit coupling 
parameter ^, while a non-cubic crystal has MCA in the 
second order in ^. Since magnetic hardness generally re- 
quires high MCA, all known hard magnets are non-cubic. 
Many of them are uniaxial, so that Aap has two principal 
components, in-plane Aab and out-of-plane Ac- 

The components of the exchange stiffness tensor may 
be found both theoretically using non-coUinear spin- 
polarized band structure calculations, and experimen- 
tally, from the long-wavelength part of the magnon dis- 
persion spectrum. Nevertheless, to my knowledge, none 
of these methods was applied to hard magnets, and the 
isotropic model {A^ — A^h — A) was explicitl y u sed in all 
micromagnetic calculations (see, e.g., Refs. I2l3l) . while A 
is usually estimated from the Curie temperature."* This 
paper reports the results of calculations of ^ab and A^ 
for several uniaxial hard magnets including CoPt, FePt, 
FePd and the SmCo5-like compound YC05. It turns out 



that exchange stiffness anisotropy in hard magnets is typ- 
ically quite large. 

Strong anisotropy of exchange stiffness may signifi- 
cantly affect the hysteretic properties of a magnet, be- 
cause it translates into anisotropic domain wall surface 
tension: 



7„=4(A„if)i/2 



(2) 



where K is the MCA constant, and An is the exchange 
stiffness along the direction normal to the domain wall. 
Therefore, the domain walls have a tendency to align 
normal to the magnetically soft direction (the one with 
the lowest An). 

On the other hand, according to the pole avoidance 
principle, the magnetostatic contribution to the free en- 
ergy prefers to eliminate the 'magnetic charges' p = 
— div M localized on the domain walls by aligning them 
parallel to the magnetization axis. For a magnet with 
easy-axis MCA, if < A^ then the exchange term 
favors the alignment of domain walls parallel to the 
easy axis, just as the magnetostatic term. However, if 
v4ab > Ac as in all studied alloys, then the exchange and 
magnetostatic terms favor different (orthogonal) domain 
wall orientations. The relative importance of these terms 
depends on the length scale and geometry of the domain 
structure. As we will see below, the exchange term often 
dominates in hard magnets with sufficiently fine domain 
structures. In this case common considerations based on 
the pole avoidance principle are inapplicable, and various 
peculiarities in the domain structure and its response to 
the external magnetic field should be expected. 

The rest of the paper is organized as follows. 
Anisotropy of exchange stiffness in several typical hard 
magnets is determined using first-principles calculations 
in Section |nj Temperature dependence of the exchange 
stiffness anisotropy is discussed in Section HTll The fol- 
lowing sections describe the effects of exchange stiff- 
ness anisotropy on the properties of the domain struc- 
ture. Anisotropic domain wall bending is considered 
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in Section IIVI from the general point of view. The ef- 
fect of anisotropic domain wall bending on the coer- 
civity of cellular Sm-Co magnets is addressed in Sec- 
tionEI Sect ion I VII examines the effects of exchange stiff- 
ness anisotropy on the structure of domain walls and co- 
ercivity of polytwinned CoPt type magnets. Section IVlII 
concludes the paper. 

II. CALCULATION OF EXCHANGE STIFFNESS 
IN CoPt TYPE AND YC05 MAGNETS 

The values of Aab, were calculated using the spin- 
spiral formulation^ of the tight-binding linear mufhn-tin 
orbiral (TB-LMTO) method within the atomic sphere 
approximation (ASA) including the 'combined correc- 
tion' term^ The spin-spiral ( "frozen magnon" ) approach 
is much more reliable compared to the calculation of ex- 
change stiffness based on a finite number of pair exchange 
parameters, because in practice the corresponding sum in 
real space does not converge. 

Local spin density approximation (LSDA) was used 
with von Barth-Hedin exchange-correlation potential^ 
The calculations were done for experimental values2i2ii£ 
of lattice constants given in Table I. The atomic sphere 
radii for both constituents were taken to be equal to each 
other in CoPt, FePt and FePd. In YC05 the radii were 
3.548 a.u. for Y and 2.627 a.u. for Co in both inequiv- 
alent positions. The calculations were repeated for two 
setups with three (^max = 2) and four (Zmax = 3) basis 
functions per atom, and care was taken to achieve con- 
vergence with the Brillouin zone sampling. Spin-orbit 
coupling was neglected. 

In the spin-spiral method the direction of the magnetic 
moment at site i depends on coordinates as 

— (sin 9 cos qr ^ , sin 9 sin qr ^ , cos 9) (3) 

where q is the wave vector and 9 the amplitude of the spin 
spiral. For a spin spiral with small q at zero temperature, 
according to JQ), we have 

E/V = Ac,f} q^qp sin^ 9. (4) 

where E and V are the excess total energy referenced 
from the ferromagnetic state and volume of the compu- 
tational cell. 

The results listed in Table ^ show that the exchange 
stiffness anisotropy is quite large in all studied alloys. 
Notably, Aab is everywhere greater than Ac due to the 
predominantly in-plane bonding between Fe or Co atoms. 
For comparison, the values of Aab and were also cal- 
culated for hep cobalt. Here the exchange stiffness is 
almost isotropic as expected, and its value of approx- 
imately 3.5x10"^ erg/cm (or D ~ 600 meV-A^) is in 
good agreement with experimenliii and with other cal- 
culations The exchange stiffness for both directions 
is larger in YC05 compared to CoPt type alloys due to 
higher Co concentration. Aab in YC05 is close to that in 
hep Co. 



TABLE I: Lattice parameters and calculated values of 
in-plane and out-of-plane exchange stiffness (units of 
10~^ erg/cm). The value given before (after) the slash was 
calculated with Zmax = 2 (Zmax = 3). 





CoPt 


FePt 


FePd 


YC05 


Co 


a, A 


3.806 


3.861 


3.860 


4.937 


2.507 


c/a 


0.968 


0.981 


0.968 


0.806 


1.623 




1.70/1.58 


1.10/0.87 


1.89/1.78 


3.97/3.89 


3.55 


Ac 


1.13/1.03 


0.38/0.06 


0.87/0.72 


1.68/1.57 


3.43 




0.66/0.65 


0.34/0.07 


0.46/0.40 


0.42/0.40 


0.97 



The out-of-plane exchange stiffness in FePt is unusu- 
ally small and very sensitive to the lattice parameter c. 
This magnetostructural effect may be quantified by the 
value W — dAc/dhic = 26 ■ 10~^ erg/cm. Low value of 
Ac and high value of W imply that moderate compres- 
sion of the order of 1-2% along the c axis may induce 
magnetic instability in FePt with the formation of a spin 
wave in the c direction. This conclusion agrees qualita- 
tively with the results of other studies suggesting that 
the layered antiferromagnetically ordered (AFM) state 
(a special case of such spin wave) in FePt has lower en- 
ergy compared to the ferromagnetic (FM) state under 
moderate c/a reductioni^ or even at experimental lattice 
parametersji^ However, strong sensitivity of to the ba- 
sis set (see Table 1) indicates that ASA is too crude for 
the description of magnetic energetics in FePt. On the 
other hand, the LSDA approximation also seems to be 
insufficient, because adding any of the two types of gra- 
dient correctionsiSiii to the LSDA exchange-correlation 
potential notably tends to stabilize the FM phase. 

Competition between different magnetic structures (in- 
cluding non-coUinear ones) is characteristic for fee phases 
of iron^'* and its alloys (FcsPt is a known Invar alloy), and 
hence it is not surprising for FePt. Indeed, while FePt is 
tetragonal, its structure fully retains the topology of the 
fee lattice if Fe and Pt sites are considered equivalent. 

From the practical point of view, structural sensitiv- 
ity of Ac in FePt suggests that the exchange stiffness 
anisotropy in this magnet may be controlled using chem- 
ical pressure, appropriate doping or off-stoichiometry. 
In view of the strong effect of this anisotropy on the hys- 
teretic properties (see below), this possibility may prove 
useful in applications, such as the design of perpendicular 
magnetic recording media. 

III. TEMPERATURE DEPENDENCE 

It is obvious that exchange stiffness anisotropy a is an 
additional parameter of micromagnetics which may have 
a strong effect on the coercivity and other properties of 
magnets. In this connection it is worth noting that in 
some materials a may strongly depend on temperature 
and doping. For example, consider a layered magnet with 
atoms of type A in even layers and type B in odd lay- 
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ers. Suppose that the exchange interaction is strong for 
A-A pairs, neghgible for B-B pairs, and small for A~B 
pairs. This is a good approximation for all CoPt type 
magnets, where A corresponds to the 3d metal, and B 
to Pt or Pd. If we assume that the magnitudes of the 
magnetic moments Ma and Mb do not depend on tem- 
perature (rigid local moments model), in the mean field 
approximation the reduced magnetization of the B layer 
TUB - {Mb) /Mb is 

ruB = f{2JBAmA/T) (5) 

where tua — (M-a) /Ma, f{x) — coth(x) — and 
Jba = X]j "^ij where site i is within the B layer and 
j runs over A sites. Since Jba, as we assumed, is much 
smaller than Jaa (defined with both i and j in the A 
layer), there is a wide range of temperatures where the 
alloy is still ferromagnetic, but Jba/T < 1. In this re- 
gion f{x) « a;/3, and, supposing that the Curie tempera- 
ture Tc is almost entirely determined by A-A interactions 
(Tc « 2JAA/i), we obtain 

ruB Jba Tc 

~ -? — , (6) 

mA Jaa J 

that is, the ratio of magnetizations of B and A sublat- 
tices is inversely proportional to temperature. Obviously, 
the relative contribution of A-B pairs to Aab and 
in the mean- field approximation follows the same law. 
At the same time, due to the layered structure the A-B 
pairs may give an important contribution to Ac at T = 
(according to the calculation, this is the case in FePt). 
Eq. © implies that close to Tc this contribution is re- 
duced by a factor Jba/ Jaa, and hence a is essentially 
determined only by exchange interaction in A-A pairs. 

In addition, in magnets like FcPt the induced magnetic 
moments of Pt atoms should be more easily destroyed by 
thermal excitations compared to the self-induced, well- 
localized Fe moments, and the temperature dependence 
of a should be even more pronounced. 

On the other hand, non-magnetic impurities in a lay- 
ered system (e.g. Cu in SmCos) may energetically prefer 
some specific layers. At high concentration of such im- 
purities the interlayer exchange coupling will be strongly 
reduced, again decreasing a and also enhancing its tem- 
perature dependence. 

These effects provide an interesting mechanism for the 
dependence of magnetic properties on temperature and 
doping due to the increased domain wall bending. This 
relationship may be important in Sm-Co type magnets 
where domain walls are heavily bent, as discussed in Sec- 
tion |V| 



IV. ANISOTROPIC DOMAIN WALL BENDING 

Bending of pinned domain walls in external mag- 
netic field was invoked by many authors to describe 



certain aspects of coercivity and hysteresis in mag- 
nets,i-— '^■^^•^^■^■^•^'^ This bending, which generally man- 
ifests itself in the initial magnetic susceptibility, was in- 
corporated in the Globus model to describe the hysteresis 
loop of granular magnets.^" Pinning of bending domain 
walls on an array of defects was also considered by a 
number of authors, see Ref. and references therein. 

In real magnets domain wall bending may play a more 
subtle role in the magnetization reversal. For some mi- 
crostructures bending of domain walls may facilitate their 
pinning by increasing the area of contact with pinning 
centers. In particular, this mechanism was discussed in 
the studies of coercivity and hysteresis loop in SmCo5 
powdersSi and in cellular Sm-Co magnetsi2i24 

In previous treatments of domain wall bending the ex- 
change stiffness was assumed to be isotropic. Here we 
will discuss the effect of exchange stiffness anisotropy on 
the domain wall bending. 

First, we will study domain wall bending neglecting 
the associated stray fields. As we will see below, this ap- 
proximation is valid when the characteristic flux closure 
length is sufficiently small. 

To get a general feeling of the problem of domain wall 
bending, it is useful to invoke a direct analogy between 
domain walls and foam bubbles. Indeed, in the pres- 
ence of external magnetic field H the domain walls in a 
uniaxial magnet experience constant pressure of magni- 
tude 2MH directed away from the regions where HM 
is positive. On the other hand, if the exchange stiffness 
tensor is isotropic, the free energy of a domain wall is 
simply proportional to its total area, just as that of a 
foam membrane. 

This analogy allows one to guess the equilibrium con- 
figurations of domain walls pinned by certain symmetric 
pinning sites. For example, a domain wall pinned by a 
circular defect should obviously have the form of a sphere 
segment with radius R related to the external field as 
R = 7(MH)~^ (the Laplace pressure 2j/Roi the curved 
domain wall compensates the applied pressure 2MH). 
This result was obtained in Ref. 24 where the sphere seg- 
ment was used as a variational trial function. A domain 
wall pinned at two parallel straight lines assumes cylin- 
drical shape with twice as smaller radius; this solution 
was discussed for a domain wall in a thin ferroelectric 
film^ 

If exchange stiffness is anisotropic, the analogy with 
the foam membrane no longer holds, because the domain 
wall surface tension is also anisotropic. Denoting the an- 
gle between the normal to the domain wall and the mag- 
netization axis as (f>, the surface tension of the domain 
wall is given by (0) with 

An = sin^ 41 + Ac cos^ 4> (7) 

Consider a domain wall pinned at a closed curve B with 
a typical size Rq lying in the x = plane, choosing the z 
axis parallel to the easy magnetization axis. The shape 
of the domain wall is defined by a function ^ = £,{y,z) 
with ^ = at the boundary B. The total free energy is 
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easily shown to be (we still neglect the stray fields and 
assume for simplicity that H is parallel to the z axis): 



F 



Tab Jl + {dyiY + - 2M7?e 



dydz (8) 



where 7ab is the surface tension of a domain wall parallel 
to the z axis, the integral is taken over the area bounded 
by _B, and dy = etc. In the isotropic case Eq. (O 

reduces to the expressions of Refs. i24M2g . In zero field 
the equilibrium shape of the domain wall is just ^ = 0, 
while at ^ it satisfies the Euler-Lagrange equation 
corresponding to the variational problem = 0: 



{dyR-^dy + ad,R-^d,) e + i?^b = 



(9) 



where R is the square root from Eq. ©, and R^^ = 
2A/i?/7ab is the domain wall curvature for "in-plane 
bending" (i.e. that with dz^ = 0). The role of exchange 
stiffness anisotropy is seen most clearly for the case when 
the external field is weak, H <^ j{MRo)~^ , and ^ -C i?o- 
In this case we obtain to first order in ^: 



{dl + adl)^ + R;,l=0 



(10) 



Eq. H10|l shows that the domain wall bends more easily 
"out of plane" {dy^ = 0) if a < 1 and "in plane" {dzS, = 
0) if a > 1. Indeed, if the boundary B is formed by 
two straight segments parallel to the z direction as in 
Ref. (e.g., pinholes in a film) at distance L = 2Rq 
from each other, then ^ does not depend on z, and from 
(|10|l we obtain the angle of domain wall deflection at the 
pinning site: /? — 2Afi?it!o/7ab- On the other hand, if 
the boundary B is formed by straight lines parallel to y 
direction (e.g., scratches on a film surface), then ^ does 
not depend on y, and the deflection angle is times 
larger. The latter conflguration is shown in Fig. ^ 

The practical implication of this result is that the ef- 
ficiency of pinning centers in a magnet with anisotropic 
exchange stiffness depends on their orientation. In a < 1 
case typical for hard magnets the domain walls bend 
more easily when pinned by defects that are normal to 
the magnetization axis. 

The exact shape of the domain wall in external field 
may be found from the solution of the Euler-Lagrange 
equation without the assumption ^ Rq. This yields 
the circular cylinder segment^ of radius i?ab for in-plane 
bending, and the elliptic cylinder segment for out-of- 
plane bending: 



(e - Co)' + a 



i?2 



ab 



(11) 



where Co = -Rab -Ro/a. 

Now, let us clarify the role of magnetostatic (stray) 
fields. If the domain wall is parallel to the magnetization 
axis z, it has no net magnetic charge. Any deviation from 
this alignment produces magnetic charge on the wall with 
surface density a = 2Mcos0, where is defined exactly 
as in Eq. Q). Obviously, in-plane domain wall bending 
does not induce any charges on the wall, and we should 



only be concerned about stray fields when we are dealing 
with out-of-plane curvature. 

To estimate when stray fields may notably affect out- 
of-plane domain wall bending, we have to compare the 
magnetostatic energy Sem generated by domain wall 
charging to the excess surface free energy Ses associated 
with this bending (both energies are defined per unit 
length in the y direction). We assume that the domain 
wall is pinned by two line defects parallel to the y axis 
and displaced from each other by a distance L = 2Ro 
along the z axis, as shown in Fig. ^ We will find Scm 
and Ses for a domain wall bent in a weak field. As 
follows from Eq. (|10|1 . in this case the domain wall is 
shaped as a segment of a circular cylinder of large ra- 
dius R — Q;i?ab- Substituting this solution in Eq. we 
obtain Seg ~ |Q;7ab^^o/-R^. 

The magnetostatic energy is given by the integral 



MHm dxdz 



(12) 



where H„j is the stray field generated by the surface 
charges on the domain wall. The distribution of these 
charges is antisymmetric with respect to the z = plane 
(see Fig.^, and the stray field obviously falls off at the 
length scale of Rq. The total positive charge per unit 
length of the wall is of the order p+ « MRq/R (assum- 
ing R ^ Rq). Since the magnetization is reversed at 
the domain wall, there is a strong cancellation in the 
integral (|12|l . Indeed, using the superposition princi- 
ple, let us add a fictitious domain wall which is a mir- 
ror image of the real domain wall with respect to the 
X = plane. This fictitious wall is shown in Fig. ^ by 
the gray dashed line. Mirror reflection also reverses the 
sign of the magnetic charges. For symmetry considera- 
tions, the contribution to the integral H12|l from outside 
of the lens-shaped area between the real and fictitious 
walls doubles when the fictitious charges are added. At 
the same time, this contribution is negligibly small, be- 
cause the two walls form a thin capacitor, and the field 
is confined to its interior. Therefore, only this interior 
region of cross-section ^Rq/R contributes to Scm- The 
stray field in this area is of the order 2p+/Ro (now we 
should take the field only from the real domain wall), and 
from lO we find Se„, « ^W^R^/R^. Thus, we obtain 
Scm/Ses = z^M^i?o/(a7ab) where the form-factor ly fa A. 
This relation also holds when the curvature is not small 
(i? ~ i?o)j but v should be somewhat different. 

The relative importance of the magnetostatic energy 
increases linearly with Rq. Using the relation 7ab — 
iKSnh/n where Jab — 7r(Aab/i^)^^^ is the in-plane do- 
main wall width, we find that 5em overcomes 5es at 
L Icr = 4a(Sab/?y, where 77 = 2ttM'^ /K is the dimension- 
less magnetostatic parameter (in hard magnets rj is small, 
e.g. in Sm2Coi7, CoPt and FePt it is close to 0.1). Thus, 
at L ^ Icr the external field works mainly against the do- 
main wall surface tension, and its anisotropy is reffected 
in the domain wall bending according to Eq. (j^. At 
L ^ ^cr the external field works mainly against magneto- 
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FIG. 1: Cylindrical out-of-plane bending of a domain wall 
(thick solid line) pinned by two line defects located at x = 0, 
z = ±Ro- C is the axis, and R the radius of the cylindrical 
domain wall. /? is the domain wall deflection angle at the 
defect. The gray dashed line shows the fictitious domain wall 
used in the calculation of the magnetostatic energy. The size 
of plus and minus symbols schematically shows the magnetic 
charge density on the real and fictitious walls. 



static forces, which makes the contribution of surface ten- 
sion (together with its anisotropy) unimportant. By min- 
imizing the total free energy given by Eq. ^ with Sem 
added, we find the deflection angle taking into account 
all energy terms: f3 ~ /3oi/(/cr + L) where Po = H/2M. 
Note, however, that this expression is approximate, be- 
cause the equilibrium shape of the domain wall at L ^ 
is no more a cylindrical segment, and because we derived 
it assuming R Rq. At L <C ?cr this result coincides 
with P found above neglecting the stray fields. However, 
as L is increased beyond Zcr, /3 approaches its asymptotic 
limit (3q and stops changing. 

In general, the parameter Zcr appears in all problems 
when magnetostatic interaction competes with the do- 
main wall energy. In a specific geometry, the character- 
istic flux closure length If should be compared with this 
parameter. At If ^ Zcr the magnetostatic interaction 
dominates; at ^/ ^ ^cr it may be neglected. In partic- 
ular, the crossover length l^ has the order of the mate- 
rial parameter A of the magnetic bubble domain theory 
(see, e.g., Ref. 0). Together with 77, this parameter also 
controls the structure of charged domain walls in thin 
filmsi^ The critical single-domain size of a particle Rsd 
is also proportional to Icr] for example, for a sphere with 
isotropic exchange stiffness i?sd ~ 5.6(5/?7>i 



V. CELLULAR Sm-Co MAGNETS 

Magnetization reversal in cellular Sm-Co magnets 
involves heavy bending of domain walls. Since the 
anisotropy of exchange stiffness strongly affects the abil- 
ity of domain walls to bend, it is likely to play an im- 
portant role in the development of coercivity. In this 
section we will focus on the model of repulsive domain 
wall pinning at the cell boundaries and show that this 
mechanism can not be realized unless the exchange stiff- 
ness is strongly anisotropic, and specifically oriented line 
defects (such as those provided by the platelet phase) are 
available for domain wall pinning in addition to the cell 
boundaries. 

Precipitation-hardened magnets based on an appropri- 
ately doped Sm2Coi7-SmCo5 system develop outstand- 
ing magnetic hardness in a wide range of temperatures. It 
is associated with the formation of a cellular microstruc- 
ture where rhomboid Sm2Coi7-based (2:17) cells are sur- 
rounded by SmCo5-based (1:5) boundary phas o^^i^^i'^° 
and is usually explained by domain wall pinning at the 
cell boundarieSi^ Some other mechanisms of coercivity 
were also suggestedi^i*^2i2^ 

Whether pinning at the cell boundaries is attractive 
or repulsive depends on the magnetic properties of 2:17 
and 1:5 phases for the given (doped) system. Recent 
experiments of Kronmiiller and Goll^ support the hy- 
pothesis that pinning is repulsive at room temperature 
but attractive at high temperatures, which also explains 
the anomalous temperature dependence of coercivity (see 
Ref. 1^ and references therein). 

The estimated unpinning field for the cell boundaries 
(in the plane-parallel configuration) agrees with the ex- 
perimentally observed coercivity. However, this does 
not fully explain high coercivity, because magnetization 
reversal always takes the path of lowest energetic barri- 
ers. It is not sufficient for a high barrier to be present; 
it is necessary that there be no way around it. The do- 
main walls might move parallel to the hexagonal axis and 
never align parallel to the cell boundaries. In order to 
be pressed against the cell boundaries, the domain walls 
must bend in the external field. SkomskiSi estimated the 
deflection angle /3 of a domain wall at a pinning site at 
56° assuming H = 0.8 T, cell size L = 80 nm, spherical 
bending, and isotropic exchange stiffness. 

However, this value of the external field appears to be 
too high for this estimate. Indeed, domain wall bending 
requires pinning at some 'seed defects' other than the cell 
boundaries/^ In order to estimate the unpinning field for 
these defects, we note that the coercivity of samples that 
had not been subjected to slow cooling is only of the order 
of 0.1 Ti^ The role of slow cooling is likely to promote 
the formation of the 1:5 phase with segregated copper^ 
Assuming that the properties of the 2:17 cells are essen- 
tially unchanged during the slow cooling, we may take 
the value of 0.1 T as an estimate of the unpinning field 
for the seed defects. The same value corresponds to do- 
main wall pinning at the vertices of the cells. Thus, we 
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obtain /3 « 6° which is clearly insufRcient to press the 
domain wall against the cell boundary. 

The above estimate assumes spherical bending and 
isotropic exchange stiffness. In reality, it is reasonable to 
assume that the values of exchange stiffness anisotropy 
in the 2:17 phases of the Sm-Co and Y-Co systems are 
very close because exchange coupling is dominated by Co- 
Co pairs. Indeed, the Curie temperatures of all R2C017 
phases (where R is a rare-earth atom or yttrium) are 
almost identical,^ while the 2:17 phase may be obtained 
from the 1:5 phase simply by a replacement of every third 
samarium atom by a C02 dumbbell. Thus, we assume 
that the factor a in pure Sm2Coi7 is close to its value 
of 0.4 obtained for YC05 at T = (Table 1). How does 
this affect the estimate of /3? As follows from Eq. Hl()(l . 
exchange stiffness anisotropy strongly facilitates domain 
wall bending only if it is pinned by line defects normal to 
the magnetization axis (as in Fig.^. In this case the de- 
flection angle /3 contains an additional factor of oT^ due 
to exchange stiffness anisotropy and a factor of 2 due to 
the fact that bending is cylindric instead of spheric. This 
brings (i to 30°, neglecting the effect of stray fields. 

The magnetostatic term in the total energy of a bent 
domain wall is notable, although not yet dominating 
for the cell size L = 80 nm. Indeed, taking Aab — 
4 X 10"^ erg/cm and a = 0.4 as found above for YC05, 
K = erg/cm^ (Ref.lH, and M = 950 emu/cm^ 

we obtain 7ab — 44 erg/cm^, i5ab = 12 nm, 7/ = 0.15, 
and Zcr ~ 130 nm. Using the approximate formula from 
Section Hvl we arrive at the final estimate of /3 « 18°. 

Let us summarize the results obtained above. For the 
typical cell size of 80 nm the external field of 0.1 T (corre- 
sponding to the unpinning of uncurved walls) may induce 
the domain wall deflection of about 18° if the exchange 
stiffness anisotropy is large and if the domain walls may 
be initially pinned by line defects normal to the magne- 
tization axis. If either of these two conditions is not met, 
the domain walls only deflect by about 6° before they are 
unpinned from the initial pinning centers. 

According to these estimates, the deflection angle does 
not reach the typical cell-boundary inclination of 30°. 
Although our assumptions may be loosened up to some 
extent (for example, allowing for a somewhat larger un- 
pinning field for uncurved walls), it seems clear that the 
repulsive cell-boundary pinning mechanism of coercivity 
may only be realized in Sm-Co magnets under a very 
favorable set of circumstances. In particular, it requires 
the presence of initial pinning defects normal to the z 
axis. 

High coercivity develops in Sm-Co magnets only when 
they are doped with zirconium.'^^ which promotes the for- 
mation of thin lamellae normal to the hexagonal axis of 
the crystal. The intersections of these lamellae with cell 
boundaries have the "right" orientation needed for high 
coercivity in the repulsive pinning case. Therefore, ij 
high coercivity of cellular Sm-Co magnets is due to re- 
pulsive pinning at the cell boundaries, it is likely that 
the lamellar phase provides the initial pinning centers 



for domain wall bending. This conclusion does not con- 
tradict the observation^ that in the low-coercivity state 
obtained after annealing at 800° C the microstructure 
and the platelet phase are fully developed. Indeed, the 
platelet phase simply provides initial pinning sites which 
may be useful only in the presence of strongly pinning 
cell boundaries developing only after the slow cooling. 
However, the microscopic origin of pinning at the initial 
pinning sites is yet to be determined. 

VI. SLANTING AND 'SMART PINNING' OF 
DOMAIN WALLS IN CoPt-TYPE MAGNETS 

In this section we will explore the effects of exchange 
stiffness anisotropy on the structure of domain walls and 
magnetization reversal in polytwinned CoPt-type mag- 
nets. The microstructure of these magnets consists of 
regular stacks of Llg-ordered domains (c-domains). In 
each stack the c-domains are separated by parallel twin 
boundaries in one of the {110} planes i^Si^Si^MMi There 
is always a high density of antiphase boundaries within 
the c-domainsMi^mia 

Usually the c-domain thickness d is large compared to 
the domain wall width 5^5 nm, and each c-domain may 
be regarded as an individual magnetic domain with in- 
trinsic 90° domain walls at the twin boundaries.'*^ The 
dynamic domain structure is formed by macrodomain 
walls^S*^ crossing many twin boundaries in a stack. 
These walls are split at the twin boundaries, and their 
segments are coupled only by relatively weak magneto- 
static forces 1^ 

Below we study the effects of exchange stiffness 
anisotropy on the properties of macrodomain walls. 
We will describe the orientation of macrodomain wall 
segments, the energetical preference of different global 
macrodomain wall orientations, the rotation of segments 
during macrodomain wall splitting in external field, and 
the relation of these properties with coercivity. 

Consider a (110) oriented macrodomain wall shown 
schematically in Fig. El As it is shown by large empty 
arrows, the magnetization in each c-domain is parallel 
to the easy c axis, and it is reversed at the macrodomain 
wall. The exchange stiffness A„ is given by Eq. |(7J where 
0, as shown in Fig. [51 is now the angle between the nor- 
mal to the domain wall segment and the tetragonal axis 
c. Without magnetostatic interaction, the surface ten- 
sion of the macrodomain wall (e.g., the free energy of the 
domain wall segment per unit normal cross-section) is 
4:{AnKy^'^ [cos{j — 0)]~*. Minimizing over </> we obtain 

tan (j) — a (13) 

where, as above, a = Ac/A^b- This result is analogous 
to a similar expression obtained for an antiphase bound- 
arytM We see that for Ac = (isotropic exchange stiff- 
ness) = 7r/4 as it should be — the domain wall has 
no preferential orientation and simply minimizes its area 
by aligning perpendicular to the twin boundaries. For 
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FIG. 2: Macrodomain wall in the polytwinned stack (three 
c-domains shown). The central segment of the wall is marked 
DW. Large empty arrows show the of magnetization. Arrows 
marked 'c' show the tetragonal axes in each c-domain; arrow 
marked 'a' shows one of the other two four-fold axes of the 
parent fee lattice. If x and y coordinates are assigned to a 
and c axes of the central c-domain, the macrodomain wall has 
a global (lIO) orientation. 

7^ ^ab the domain wall segments slant so as to de- 
crease the surface tension. 

Segments of a (110) macrodomain wall carry magnetic 
charges of alternating signs, "'"^ and the flux closure length 
I f is obviously of the order of the c-domain thickness d. 
Therefore, at d ^ the magnetostatic interaction dom- 
inates, and the segments align parallel to the tetragonal 
axis to get rid of the magnetic charge, while at d ^ 
their orientation is determined by Eq. H13|) . Here Icr is 
proportional to S/rj, but the form-factor is different from 
that found in Section IIVI for domain wall bending. In 
CoPt and FePt 4r is of the order of 50 nm. Thus, while 
d increases during the microstructural coarsening, the do- 
main wall segments gradually slant from the angle given 
by towards the tetragonal axis.^^ 

Contrary to the (110) oriented macrodomain wall, the 
segments of a (001) oriented one are perpendicular to 
the twin boundaries at any d, because this minimizes 
both the surface energy and the magnetostatic energy. 
At d > l„ the (001) macrodomain wall has lower energy 
than the (110) one because its segments do not carry 
any magnetic charge. However, at c? ^ Icr the prefer- 
ential global orientation is determined by exchange stiff- 
ness anisotropy. The surface tension for (001) and (110) 
macrodomain walls (per unit normal cross-section) is 7ab 
and 7ab [2a/{l + a)]^/^, respectively. Therefore, at a > 1 
the (001) orientation is favorable, while in the typical 
case a < 1 the (110) orientation is favorable. This ex- 
plains why macrodomain walls observed in FePt crystals 
have {110} orientation^ while those in FePd with much 
larger rj and hence smaller l^r have {100} orientation."'^ 

High coercivity of CoPt-type magnets is likely due to 
the combination of macrodomain wall splitting and pin- 
ning of their segments at antiphase boundaries4£i^ The 
highest possible coercivity is achieved when antiphase 
boundaries are planar, and domain wall segments are par- 
allel to them."*- The antiphase boundaries in CoPt-type 



magnets often have a preferential crystallographic orien- 
tation, which is clear from theoretical considerations,^^ 
and observed experimentally for CoPti^ Typically they 
slant towards the tetragonal axis, contrary to the domain 
wall segments which slant away from the tetragonal axis 
at Ac < Aah and d <g; Icr- This means that while d in- 
creases during the microstructural coarsening, at some 
point the domain wall segments will become parallel to 
the preferential orientation of antiphase boundaries pro- 
ducing a maximum in the coercivity. 

Interestingly, a similar segment rotation may occur 
dynamically. If the external field above the splitting 
threshold^ is applied parallel to the twin boundaries, the 
macrodomain wall is split in two "partial macrodomain 
walls" which are driven apart from each other. One par- 
tial wall is composed of segments in all odd c-domains, 
and the other of segments in even c-domains. Suppose 
that d <^ Ici so that the orientation of segments in the 
(110) macrodomain wall is given by Eq. (|13() . Since 
these segments carry magnetic charges of ±2M cos per 
unit area, the two partial macrodomain walls carry equal 
charge densities of opposite sign. The corresponding 
stray field makes an additional contribution AE to the 
magnetostatic energy proportional to the distance L be- 
tween the partial macrodomain walls. If the angle 4> were 
fixed, at L > /cr this positive contribution would dom- 
inate over the surface energy of the segments (now the 
flux closure length If is clearly L). Therefore, as the two 
partial macrodomain walls move apart (L is increased to 
Icr), their segments gradually rotate toward the tetrago- 
nal axis to get rid of the magnetic charge. 

As a result, at some L the segments of a splitting 
macrodomain wall become parallel to antiphase bound- 
aries, just as in the case of increasing d discussed above. 
If the typical distance between antiphase boundaries is 
smaller than Icr (the scale of L where domain wall seg- 
ment rotation occurs), the segments will be pinned by 
antiphase boundaries at the plain-parallel configuration. 
In this scenario the coercivity achieves its highest possi- 
ble value for suitably oriented polytwinned stacks at any 
d Icr, i-e- at relatively early stages of coarsening. This 
mechanism may play an important role in real CoPt-type 
magnets developing high coercivity just at these early 
stages. 



VII. CONCLUSION 

Using the spin-spiral version of the TB-LMTO method, 
the in-plane and out-of-plane principal components of the 
exchange stiffness tensor were calculated for several typi- 
cal hard magnets. The results show that this tensor usu- 
ally has a considerable anisotropy. The out-of-plane com- 
ponent is smaller than the in-plane one in all studied hard 
magnets except pure hep Co where exchange stiffness is 
isotropic. The anisotropy is especially high in FePt. In 
certain materials with intrinsically non-magnetic layers 
(CoPt, FePt, etc.) the anisotropy of exchange stiffness 
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may strongly increase at finite temperatures or with suit- 
able non-magnetic doping. 

Anisotropy of exchange stiffness may have a strong ef- 
fect on the orientation of domain walls and on their resis- 
tance to bending, and hence on the hysteretic properties 
of the magnet. These effects may be expected whenever 
the typical flux closure length associated with the stray 
fields does not exceed the crossover length Icr (x S/r] (see 
Section Hvjl. 

Low out-of-plane exchange stiffness facilitates out-of- 
plane domain wall bending. This effect is crucial for the 
development of coercivity in cellular Sm-Co magnets in 
the repulsive cell-boundary pinning regime. This regime 
may be realized only in the presence of linear pinning 
defects normal to the magnetization axis (such as the 
intersections of lamellae with the cell boundaries) pro- 
viding initial pinning centers necessary for the domain 
walls to bend and get pressed against the strongly pin- 
ning cell boundaries. However, the estimates obtained 
in Section suggest that the set of conditions for the 
realization of this regime is very strict. 

In polytwinned CoPt-type magnets the competition 
between the anisotropy of exchange stiffness and mag- 
netostatic interaction controls the orientation of domain 
wall segments and the preferential global macrodomain 
wall orientation. In particular, the domain wall seg- 
ments gradually rotate toward the tetragonal axis during 
the microstructural coarsening as the c-domains become 



thicker. If the antiphase boundaries have a preferential 
orientation, the coercivity achieves its maximum at the c- 
domain thickness when the domain wall segments become 
parallel to the antiphase boundaries. The same competi- 
tion also leads to the dynamic rotation of the segments of 
a macrodomain wall which is being split by the external 
field. This rotation may result in a dynamic self-locking 
of domain wall segments at antiphase boundaries. 

The effects discussed in this paper demonstrate that 
the competition between exchange stiffness anisotropy 
and magnetostatic energy is a crucial driving force be- 
hind the formation and dynamics of nanoscale domain 
structures in hard magnets. 
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